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ABSTRACT 


In  this  thesis,  we  emulate  a  Wind  Turbine  Generator  by  driving  a  Doubly  Fed  Induction 
Generator  (DFIG)  via  a  DC  motor  with  variable  input  torque  capability.  The  two  circuits 
of  concern  are  the  DFIG  and  Supply-side  circuits.  They  are  electrically  coupled  with 
back-to-back  Space  Vector  Modulators  (SVMs)  that  are  coupled  through  a  DC  Link 
consisting  of  a  capacitor  bank.  The  SVMs,  with  DC  Link,  provide  bi-directional  power 
flow  between  the  DFIG  rotor  and  Supply-side  power  supply. 

The  Supply-side  circuit  senses  voltage,  current  and  DC  Link  voltage  (VDC),  and 
sends  the  sensed  inputs  to  a  Field  Programmable  Gate  Array  (FPGA),  which  is  used  to 
derive  control  of  VDC  so  that  it  is  maintained  at  200V  via  SVM. 

The  DFIG  circuit  senses  rotor  current,  stator  voltage,  rotor  speed,  and  rotor 
electrical  position,  which  are  used  to  derive  a  current  reference  that  is  used  to  control  the 
rotor  speed  using  a  Proportional  Integral  Controller  (PI).  For  our  application  we  show 
that  the  DFIG  rotor  speed  can  be  controlled  (subsynchronous,  synchronous,  or 
supersynchronous)  for  a  given  input  torque  and  thus  distributes  the  power  generated  by 
the  DFIG  as  predicted  between  the  rotor  and  stator. 
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EXECUTIVE  SUMMARY 


In  this  thesis,  we  emulate  a  Wind  Turbine  Generator  by  driving  a  Doubly  Fed 
Induetion  Generator  (DFIG)  via  a  DC  motor  with  variable  input  torque  eapability. 
Previous  researeh  has  shown  we  ean  use  a  DFIG  with  baek-to-back  SVMs  to  aeeomplish 
a  deeoupled  Supply-side,  and  DFIG-side,  eontrol  seheme  while  allowing  power  flow  to 
occur  between  the  two  systems. 

The  two  circuits  of  concern  are  the  DFIG  and  Supply-side  circuits.  The  DFIG  and 
Supply-side  circuits  are  electrically  coupled  with  back-to-back  (DFIG-side  and  Supply- 
side)  Space  Vector  Modulators  (SVMs)  that  are  coupled  through  a  DC  Link  consisting  of 
a  capacitor  bank.  The  back-to-back  SVMs,  with  DC  Link,  provide  bi-directional  power 
flow  between  the  DFIG  rotor  and  Supply-side  power  supply.  Bi-directional  power  flow  is 
achieved  when  the  DFIG  is  controlled  at  supersynchronous  speeds,  and  input  torque  is 
enough  to  overcome  losses,  and  subsynchronous  speeds  (Supply-side  provides  power  to 
the  rotor). 

The  Supply-side  power  supply  (3-phase  AC)  provides  power  to  the  Supply-side 
circuit  (stepped  down  by  a  VARIAC)  and  the  Stator  of  the  DFIG.  The  Supply-side  circuit 
senses  the  supply  side  voltage,  current  and  DC  Link  voltage  (VDC)  and  sends  the  sensed 
inputs  to  a  Field  Programmable  Gate  Array  (FPGA).  The  FPGA  is  programmed  via  a 
Simulation  software  package  called  Simulink  to  achieve  the  desired  control  of  the 
specified  parameters.  For  the  Supply-side  circuit  we  use  the  mentioned  sensed  inputs  to 
derive  control  of  VDC  so  that  it  is  maintained  at  200V  via  SVM.  There  is  also  a  built-in 
protection  feature  in  the  Supply-side  circuit  that  will  turn  off  the  SVM  if  a  VDC  of  240V 
is  reached. 

The  DFIG  circuit  power  supply  is  provided  by  the  DC  Link  and  DFIG-side  SVM. 
The  DFIG  circuit  senses  rotor  current,  stator  voltage,  rotor  speed,  and  rotor  electrical 
position  (as  provided  by  an  encoder).  These  inputs,  like  the  Supply-side  circuit,  are  then 
sent  to  an  FPGA  which  is  programmed  to  control  the  DFIG  parameters.  The  DFIG  circuit 
will  take  the  rotor  speed  input  and  derive  a  current  reference  that  is  used  to  control  the 
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rotor  speed  using  a  Proportional  Integral  Controller  (PI).  It  is  important  to  show  that  we 
can  control  the  DFIG  rotor  speed  since  traditional  Wind  Turbine  Generators  would  use 
turbine  blade  pitch  control  in  tandem  with  a  synchronous  generator  to  provide  a  constant 
speed,  and  hence  a  constant  output  frequency.  Further,  it  is  noteworthy  that  a  wind-speed 
profile  can  be  adapted  to  a  specific  Wind  Turbine  (outfitted  with  a  DFIG)  in  order  to 
program  the  DFIG  to  run  at  specified  speeds  for  a  given  wind  speed.  This  tool  can  be 
used  to  provide  an  optimized  rotor  speed  for  a  given  wind  speed.  For  our  application  we 
will  show  that  the  DFIG  rotor  speed  can  be  controlled  (subsynchronous,  synchronous  or 
supersynchronous)  for  any  given  input  torque  and  thus  distributes  the  power  generated  by 
the  DFIG  as  predicted  between  the  rotor  and  stator.  Further,  this  experiment  will 
demonstrate  the  bi-directional  power  flow  of  the  back-to-back  SVMs  showing  power 
recovery  via  the  DFIG  rotor. 

As  illustrated  in  the  below  figure,  the  Supply-side  SVM  converts  a  3 -phase  AC 
source  to  a  constant  DC  Link  voltage  of  200V.  The  DC  Link  voltage  is  then  inverted  into 
a  3-phase  waveform  via  the  DFIG-Side  SVM  (during  modes  for  which  power  is  being 
delivered  to  the  DFIG).  The  digital  control  for  the  Insulated  Gate  Bipolar  Transistor 
(IGBT)  network  is  implemented  via  the  FPGA.  The  FPGA  interface  samples  the  rotor 
current  and  the  stator  voltage  of  the  DFIG  through  an  analog  to  digital  (A/D)  converter. 
There  is  an  input  from  the  encoder,  Or,  to  provide  the  position  and  speed  of  the  rotor  to 
the  FPGA  through  the  interface  circuit  board.  An  algorithm  then  calculates  the  slip 
frequency  of  the  DFIG  which  is  used  to  derive  the  needed  rotor  current  frequency  and 
amplitude  to  maintain  a  specified  rotor  speed,  and  therefore  a  controllable  output 
frequency. 
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Simplified  Circuit  of  Wind  Turbine  DFIG  System. 


The  below  figure  demonstrates  an  input  torque  transient  representing  a  large  wind 
gust.  By  commanding  a  specified  supersynchronous  speed,  while  taking  into  account 
windage  and  friction  losses  of  the  system,  it  is  clear  that  the  input  torque  transient 
produces  a  power  flow  reversal  in  the  rotor.  In  other  words,  the  rotor  mode  transitions 
from  power  flow  from  the  rotor  (negative  rotor  current)  to  power  flow  to  the  rotor 
(positive  rotor  current). 

In  addition,  this  transient  plot  (simulation  model)  shows  a  smooth  transition  as  the 
power  flow  reverses  between  power  flow  from  the  rotor  to  power  flow  to  the  rotor.  This 
transition  is  modeled  using  a  constant  commanded  rotor  speed  along  with  a  10%  decrease 
in  input  torque.  The  next  two  figures  (experimental  results)  parallels  this  transition  with 
the  steady  state  captures  corresponding  to  the  simulation  transient.  The  experimental 
results  were  run  through  a  low  pass  filter  with  a  comer  frequency  of  80  Hz  in  order  to 
create  a  useful  plot  of  the  event.  Due  to  windage  and  friction  losses,  a  commanded 
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supersynchronous  speed  was  ordered  to  achieve  the  power  flow  transition  results  during 
the  constant  rotor  speed  experiment.  If  the  experimental  results  are  closely  examined,  it  is 
interesting  to  note  the  small,  low  frequency  ripple  due  to  speed  control  beating. 
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Simulation:  Input  Torque  Transient  Showing  Rotor  Power  Reversal. 
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Experimental:  Rotor  Current  and  Voltage  with  Power  Flow  from  the  Rotor. 
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In  conclusion,  we  showed  that  we  can  control  the  DFIG  rotor  excitation  to 
provide  a  constant  output  frequency.  Our  application  demonstrated  that  the  DFIG  rotor 
speed  can  be  controlled  (subsynchronous,  synchronous  or  supersynchronous)  for  a  given 
input  torque  and  thus  distributes  the  power  generated  by  the  DFIG  as  predicted  between 
the  rotor  and  stator.  Finally,  this  experiment  also  demonstrates  the  bi-directional  power 
flow  of  the  back-to-back  SVMs  showing  power  recovery  via  the  DFIG  rotor. 
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I. 


INTRODUCTION 


A.  MISSION 

In  this  thesis,  we  emulate  a  Wind  Turbine  Generator  by  driving  a  Doubly  Fed 
Induction  Generator  (DFIG)  via  a  DC  motor  with  variable  input  torque  capability.  As 
discussed  in  [1-4]  we  can  use  a  DFIG  with  back-to-back  SVMs  to  accomplish  a 
decoupled  Supply-side,  and  DFIG-side,  control  scheme  while  allowing  power  flow  to 
occur  between  the  two  systems.  The  control  schemes  developed  with  this  research  can 
branch  into  shipboard  application  as  discussed  below. 

The  benefits  of  researching  and  developing  a  Wind  Turbine  Generation  System, 
refitted  with  a  DFIG,  span  a  wide  range  of  purposes  for  the  United  States  as  a  whole,  as 
well  as  the  United  States  military.  In  particular,  there  are  several  naval  warships  still 
running  large,  separately  supplied,  gas  turbine  generators  to  provide  the  ship  with 
electrical  power  to  vital  loads.  Considering  this  one  example,  a  properly  researched  and 
developed  DFIG  with  precision  characteristic  control  can  not  only  save  space  (in  a  space¬ 
scarce  environment),  but  also  provide  a  more  efficient  means  of  power  distribution  when 
mechanically  coupled  to  the  main  propulsion  gas  turbine  shaft  line.  This,  of  course, 
requires  a  variable  speed  constant  frequency  control  scheme  since  the  propulsion  turbine 
shaft  speed  is  variable  according  to  the  desired  ship  speed  [5].  The  applications  of  the 
DFIG  are  many  and  complex,  but  the  simple  example  above  shows  a  common  sense 
approach  to  save  space  and  fuel  consumption  for  our  United  States  military. 

B.  OBJECTIVE 

The  goal  of  this  thesis  was  to  build  a  computer  simulation  of  a  Wind  Turbine 
Generation  System,  based  on  mathematical  equations,  through  the  use  of  Mathworks’ 
Simulink  software.  Once  the  simulation  was  complete,  the  goal  shifted  to  build  a  fully 
functioning  model  of  a  Wind  Turbine  Generation  System.  A  series  of  experiments  were 
performed  to  compare  actual  results  with  expected/simulation  results.  Once  the  series  of 
experiments  proved  the  simulation  was  an  accurate  depiction  of  the  model,  the  machinery 
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was  put  to  full  test  as  a  Wind  Turbine  Generation  System  to  be  used  in  future  lab  work 
and  research  to  hone  in  features  that  could  be  discovered. 

C.  APPROACH 

The  first  step  is  to  generate  a  computer  simulation  of  the  system,  shown  in  Figure 
1.  There  are  a  number  of  ways  to  simulate  the  DFIG  coupled  with  a  Supply-side,  via 
back-to-back  SVMs,  but  the  one  used  for  this  thesis  is  covered  in  [1],  Specifically,  a 
computer  model  of  the  DFlG-side  circuit  based  on  mathematical  equations  is  simulated  in 
the  rotor  reference  frame  using  dqO  transformation.  The  advantage  of  using  dqO 
transformation  is  most  apparent  when  attempting  to  decouple  DFlGs  real  and  reactive 
power.  By  doing  so,  it  is  possible  to  separately  control  generator  torque  (hence  speed) 
and  reactive  power  (hence  rotor  and  stator  power  factor)  when  the  currents  controlled  are 
in  the  slip  frequency  reference  frame.  The  DFlG-side  circuit  simulation  inputs  are  input 
torque,  stator  voltage,  rotor  current,  and  rotor  speed  and  position. 
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Figure  1 .  Simplified  Circuit  of  Wind  Turbine  DFIG  System. 

The  Supply-side  circuit  simulation  is  based  in  the  stator  reference  frame  using  dqO 
transformation.  Not  unlike  the  DFIG-side  circuit,  there  is  a  distinct  advantage  in  using  the 
dqO  transformation  so  as  to  decouple  real  and  reactive  power  on  the  Supply-side  circuit. 
This,  in  turn,  allows  a  separate  control  of  DC  Link  voltage  (via  a  current  reference)  and 
reactive  power  (hence  Supply-side  power  factor).  The  Supply-side  circuit  simulation 
inputs  are  Supply-side  voltage,  DC  Link  voltage  and  Supply-side  current. 

The  physical  system  consists  of  a  DFIG,  back-to-back  SVMs,  Student  Design 
Center  (SDC)  [6],  an  encoder,  and  a  DC  motor  used  to  simulate  changes  in  wind  speed. 
The  DFIG  is  connected  to  the  grid  via  the  stator  windings  and  to  the  DFIG-side  SVM  via 
the  rotor  windings.  The  SDC  is  the  mechanism  by  which  the  control  algorithm  is 
implemented.  The  SDC  measures  the  applicable  system  parameters  and  outputs  the 
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appropriate  control  signals  to  the  SVM  to  control  the  rotor  currents  in  a  fashion  that 
allows  for  a  specific  mode  of  operation.  For  our  application  we  will  show  that  the  DFIG 
rotor  speed  can  be  controlled  (subsynchronous,  synchronous  or  supersynchronous)  for  a 
given  input  torque  and  thus  distributes  the  power  generated  by  the  DFIG  as  predicted 
between  the  rotor  and  stator. 

D.  THESIS  ORGANIZATION 

This  thesis  is  organized  with  Chapter  II  covering  the  theory  of  induction 
machines;  to  include  the  equations  that  represent  voltage,  flux  linkage,  and  reference 
frame  transformations  and  how  the  voltage  and  flux  equations  and  machine  parameters 
are  used  to  model  the  system.  Chapter  III  is  dedicated  to  the  control  algorithm.  It  explains 
how  an  algorithm  is  formulated,  based  on  the  equations  derived  in  Chapter  II,  and  is 
manipulated  to  control  the  simulation  and  the  physical  system.  Chapter  IV  displays  the 
results  of  the  working  system,  and  compares  the  actual  results  to  the  simulation  results. 
Chapter  IV  also  includes  the  testing  procedure  used  to  calibrate  and  implement  the 
encoder  and  its  offset,  the  method  used  to  transform  the  DFIG-side  rotor  currents  to  a  slip 
frequency  reference  frame,  and  Xilinx  implementation  methods.  Chapter  V  discusses  the 
conclusions  and  future  research  objectives  for  this  thesis.  The  Appendix  contains  Matlab 
files,  data  sheets  for  equipment  used,  specific  schematic  layout  of  the  Simulink  systems, 
and  transformation  derivations  used  in  this  thesis. 
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II.  DFIG  THEORY 


A.  INDUCTION  MACHINE  EQUATIONS 

Simulink  is  a  powerful  tool  that  can  be  implemented  in  the  simulation  of  any 
given  system  that  can  be  mathematically  represented.  To  understand  the  simulation  of  the 
Wind  Turbine  Generation  System,  refitted  with  a  DFIG,  there  must  be  a  thorough 
understanding  of  the  equations  that  represent  the  DFIG.  The  equations  used  to 
characterize  the  DFIG  are  obtained  from  [7].  The  induction  machine  voltage  equations 
are  expressed  as 

V’abcs]  =  Vs][hbcsV  p[Kbcs] 

V’abcr]=[rr]VabcrV  P[\bcr] 

where  the  subscript  s  and  r  are  denoting  the  stator  and  rotor  associated  variables  and  X 
represents  the  flux  linkage.  Finally,  p  is  used  as  the  derivative  operator.  For  a 
magnetically  linear  system,  flux  linkage  can  be  represented  as  equation  (2). 
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where  Lsr  represents  the  mutual  induetanee  between  the  stator  and  rotor.  The  stator,  rotor 
and  mutual  winding  induetanees  are 
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(5) 


where  the  variables  and  Lms  are  the  leakage  and  magnetizing  inductances  of  the  stator 
windings;  conversely,  Lir  and  L^r  are  the  leakage  and  magnetizing  inductances  for  the 
rotor  windings.  Lastly,  0r  is  the  angular  position  of  the  rotor. 


When  creating  an  equivalent  circuit,  a  standard  convention  is  to  refer  the  rotor 
variables  to  the  stator  windings  using  the  appropriate  turn  ratio  N. 


cibcr  ducr 


^abcr  ^abcr 


^  N 


(6) 

(7) 

(8) 

(9) 


This  turns  ratio  referral  also  applies  to  the  magnetizing  and  mutual  inductances  and  are 
related  by 


L  =^L 

ms  SI 


L'  = 


K  = 


Jlrj 


(10) 


(11) 


(12) 


N 

where  by  defining  Ll  =^4^and 


^  N 


the  rotor  inductance  variables  can  be 


referred  to  the  stator  windings  as  shown  in  equations  (13)  and  (14). 
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(14) 


The  flux  linkage  and  voltage  equations  in  (1)  and  (2),  in  terms  of  variables  referred  to  the 
stator  windings,  ean  now  be  expressed  as 


^abcs 

1 - 

1 _ 

^abcs 

2  ' 

_  ^abcr  _ 

•  f 

J'abcr  _ 

V’abcs]  =  Vs]Vabcs]  +  P[Kbcs] 
V’'abcr]  =  [K]VabcrV  p[Kbcr] 


(15) 


(16) 


Above,  the  physieal  voltage  equations  were  given  for  the  induetion  machine;  however,  in 
an  effort  to  reduce  the  complexity  of  their  differential  equations  we  will  use  a  change  of 
variables  transformation  to  a  new  reference  frame  as  described  in  equation  (17). 


(17) 


where  the  matrix  for  stator  variables  and  for  rotor  variables  is  expressed  as 
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cos  (9  cos(^-^)  cos(^  +  ^) 

K^=—  sin^  sin(^  +  ^) 

ill 

_  2  2  2  . 

(18) 

COSy5  COs(/?-^)  COs(y5  +  ^) 

K^=—  sin  J3  sin  (y^  -  sin  (yff  + 

11  1 

.  2  2  2  . 

p  =  e-0^ 

(19) 

In  equations  (18)  and  (19),  ^  is  the  variable  used  to  transform  the  variables  for 

which  equation  (17)  are  referred.  The  place  holder  x  in  equation  (17)  represents  this 
particular  frame  of  reference  that  variables  are  transformed.  The  place  holder  / 
represents  the  voltage,  current,  or  flux  linkages  in  the  stator  or  rotor  windings.  If  no  place 
holder  is  used  then  the  reference  frame  is  considered  to  be  an  arbitrary  frame  of 
reference.  When  transforming  between  reference  frames  equation  (20)  is  utilized  where 
y  represents  the  new  frame  of  reference. 

cos(6»^-^J  sm[e^-e^)  0 

cos(^^-^4  0 
0  0  1 

(20) 

After  performing  the  transformation  of  equation  (17)  on  each  component  in  the  matrices 
on  the  right  hand  side  of  equation  (16),  the  voltage  equations  can  now  be  expressed  in  the 
arbitrary  reference  frame  as 
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\s=rshs+PKr 

^',r=r!i[r+{^-^r)Kr^PKr 

^'dr=r'/jr-{0)-(0,)K,+pKr 

^Or  ~  ^rhr  “*“  P\r 

(21) 


In  equation  (21),  <y^is  electrical  angular  frequency  of  the  rotor,  and  CO  is  the  angular 
frequency  of  the  reference  frame  that  the  variables  are  transformed  in  radians  per  second. 
Since  y/  =  and  co^^  =  2;7-60  ™54c  ^  it  is  convenient  to  express  the  inductance  component 

of  the  flux  linkages  in  equation  (21)  as  a  reactance  instead  of  as  an  inductance,  beginning 
with  equation  (22). 


0)  p 

a,  CO, 

CO  p 

=%* - ¥qs+—¥ds 

CO,  CO, 


a, 

,  ,  P  , 


qr  r  qr 


CO,  CO, 


qr 


^dr  =  - 


COu 


COu 


^Or  ^rhr  ¥or 

co^ 


(22) 


The  flux  linkage  component  of  equation  (21)  now  has  units  of  flux  linkage  per  second 
and  represented  by  the  symbol^ .  For  ease  of  understanding  the  flux  linkages  in  (21)  and 
flux  linkage  per  second  in  (22)  are  compared  in  equations  (23)  and  (24). 
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—  Li  i  L  (i  i  ^ 

qs  Is  qs  ms\  qs  qr  ^ 

^ds  ~  ^Ishs  ^ms  ^hs  hr  ) 

\s  ~  ^Ishs 

A'  =  Lj  i  +L  ii  +/'  ) 

Kr  =  Khr  +  ^ms  {L  +  4  ) 

0  '  _  T  j'f 

^Or  ~  Hrhr 

(23) 

¥qs  =  ^Jqs  +  ^ms  ^  C  ^ 

¥ds  =  ^IsL  +  (4  +  ^'dr  ) 

¥os  = 

¥'qr  =  ^'idqr  +  ^ ms  V  ) 

y^'dr  =  Kddr  +  ^ms  (4  +  4  ) 

¥l=X'Ar 

(24) 

The  above  equations  mathematically  represent  the  behavior  of  a  DFIG,  and  are  used  in 
the  simulation  process  to  arrive  at  an  appropriate  physical  model  that  will  be  discussed 
below.  More  can  be  read  concerning  induction  machine  theory  as  used  in  this  thesis 
application  in  [4], 

B.  SYSTEM  CONFIGURATION 

The  combined  system  is  shown  in  Figure  1.  The  DFIG  is  a  175W,  120V,  60Hz,  4- 
pole  machine  with  parameters  listed  in  the  Appendix  data  sheet.  The  system  parameters 
were  used  in  order  to  effectively  model  the  system  in  Simulink.  The  Lab-Volt  model 
823 1  was  used  as  the  DFIG.  The  stator  consists  of  3-phase  wye-connected  windings  with 
a  turns  ratio  of  Ns/Nr  =  10  and  resistance  Rs  =  I2Q.  The  rotor  consists  of  3-phase  wye- 
connected  windings  with  resistance  Ry  =  4Q.  The  reactance  of  the  DFIG  are  Xm=  I80Q 
dindXis=Xir=  90.. 

The  DC  motor  used  to  emulate  the  wind  turbine  is  the  Lab- Volt  8211.  It  is  a 
I75W,  I20V,  4-pole  machine  with  a  nominal  operating  speed  of  1800  RPM.  The 
operating  speed  range  used  was  1350  RPM  to  2150  RPM. 
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The  power  supply  used  for  the  experiments  was  the  Lab-Volt  8821-10.  It  is  a  3- 
phase  120V  to  208V  AC  power  supply  with  a  separate  VARIAC  power  supply  and  120V 
DC  power  supply.  In  order  to  properly  provide  power  to  the  Supply-side  circuit,  the 
VARIAC  output  is  used  for  the  Supply-side  SVM  with  a  nominal  setting  of  60V  rms. 
This  allows  sufficient  voltage  to  be  boosted  to  maintain  the  DC  Link  at  200V  without 
excessive  current  transients.  Further,  a  set  of  640pH  inductors  were  placed  in  series  (see 
Figure  2.)  to  the  input  of  the  SVM  in  order  to  decouple  the  current  sensors  in  the  SDC 
from  the  switching  action  of  the  SVM. 

l-r  Rr 


Figure  2.  Simplified  Single  Phase  Circuit  Showing  the  Inductor  Connections. 


The  mathematical  analysis  used  to  arrive  at  the  equivalent  simulation  matrices 
used  in  Simulink  are  shown  in  equations  (25)-(31).  The  single  phase  neutral  point 
voltages  are 


(25) 


and  the  KVL  analysis  results  in 


Kn2  +  -  Kanl  +  (^1  ~ 


(26) 
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Knl  +  +  Rfh  -  Kbnl  +  (^1  ~  ^2) 


(27) 


Knl  +  +  R/c  -  Kcnl  +  (^1  ~  ^2) 


(28) 


Where  nj  and  n2  are  neutral  point  voltages. 

Subtracting  (26)  and  (27)  gives 

^an2  ~  ^bn2  (4  ~  4  )  “*"  (4  ~  4  )  ~  ^sanl  ~ 

(29) 


and  subtracting  (27)  and  (28)  gives 


^bn2  ^cn2'^^^fih  4)'*''^/'(4  h^~^sbnl  K 


scnl 


(30) 


Defining  ,  and  4  +  i,  +4=0  the  matrix  yields 
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Next,  equation  (31)  can  be  put  into  a  state  space  form  by  defining 
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So  that  (31)  can  be  written  as 

Mx  +  sNx  +  Pu  =  0. 

(32) 

Then  the  state  space  model  becomes 

+  Bu 

(33) 

with 


A  =  -N~'M 
B  =  -N~'P 

(34) 
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Therefore,  the  matrices  used  in  Simulink  are 
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(35) 

The  above  process  represents  the  common  approach  used  to  parallel  circuit 
analysis  and  Simulink  modeling.  Since  Simulink  is  a  mathematical  based  simulation 
program,  vice  any  given  circuit  solving  based  program,  it  is  imperative  to  accurately 
analyze  the  mathematical  characteristics  of  any  given  subsystem  in  order  to  effectively 
implement  Simulink. 

Once  an  in  depth  understanding  of  induction  machine  theory  and  system 
configuration  was  established,  specific  control  methods  were  explored.  The  control 
methods  included  control  components,  Supply-side  control,  DFIG-side  control,  SVM 
scheme  concepts,  and  the  encoder  implementation. 
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III.  CONTROLLER  SCHEME 


A.  CONTROLLER  COMPONENTS 

The  back-to-back  SVM  components  used  for  this  research  are  a  combined 
package  manufactured  by  SEMIKRON.  The  package  is  called  the  SEMISTACK-IGBT 
and  consists  of  a  three  phase  inverter  that  is  coupled  to  the  DC  Link  capacitor  bank.  The 
inverters  inside  the  SEMISTACK  are  made  of  SKM  50  GB  123D  IGBTs  that  are 
controlled  by  SEMIKRON  SKHI  22  gate  drivers.  The  Supply-side  SEMISTACK  is  used 
to  maintain  the  DC  Link  at  200V  and  can  also  be  used  to  control  the  power  factor  to  the 
grid  when  in  inverter  mode.  Controlling  the  grid  power  factor  was  not  within  the  scope  of 
this  thesis,  outside  of  experimentally  showing  its  explicit  control,  but  is  a  feature  that  can 
be  used  to  minimize  transmission  losses  in  a  large  scale  project.  The  control  for  each 
IGBT  in  the  SEMISTACK  is  implemented  through  the  SDC.  Two  SDCs  are  utilized 
(Supply-side  and  DFIG-side)  and  include  a  Xilinx  FPGA  and  a  Control  board  that  has 
various  inputs  and  outputs  that  interface  with  the  Supply-side  and  DFIG-side  circuits. 
One  of  the  inputs  is  an  8-channel,  12-bit  A/D  converter  used  to  measure  the  voltage  and 
current  signals  from  their  associated  systems.  The  DFIG-side  circuit  also  has  a  4-pin 
input  that  samples  the  rotor  encoder  pulse  waveforms  and  converts  them  into  a  rotor 
position  and  rotor  speed  input  via  an  algorithm  implemented  in  SIMULINK.  The  FPGA 
is  programmed  to  take  the  sampled  inputs  and  generate  the  desired  gate  signals  to  the 
control  board.  The  control  board  has  BNC  (Bayonette  Neil-Concealman)  connecters  that 
connect  to  the  SEMISTACK’ s  gate  drivers.  The  gate  drive  signals  activate  the  IGBTs  in 
a  fashion  that  will  produce  the  desired  3 -phase  current  in  the  rotor  windings  for  the 
DFIG-side,  or  the  commanded  DC  Link  voltage  for  the  Supply-side.  The  procedure  for 
using  the  SDC,  architecture  of  internal  components,  and  the  benefits  of  using  FPGA  for 
digital  control  of  power  systems  is  discussed  in  [6]. 

B.  SUPPLY-SIDE  CONTROLLER 

The  Supply-side  circuit  is  powered  via  a  3-phase  AC  power  supply  stepped  down 

by  a  VARIAC  to  60V  rms.  Of  note,  it  is  important  to  calculate  the  DC  Link  injection 
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voltage  as  described  in  [2],  The  injected  voltage  should  be  less  than  /  2V2  , 

which  corresponds  to  an  injected  voltage  of  60V  rms.  There  is  a  set  of  three  inductors 
connected  in  series  to  the  SEMISTACK,  downstream  of  the  Supply-side  SDC,  which 
decouples  the  SDC  sensors  from  the  switching  action  of  the  IGBTs.  Figure  1.  shows,  and 
equations  (25)  through  (35)  describe  their  matrix  analysis.  The  SDC  senses  the  stepped 
down  power  supply  line-to-line  voltages,  Vab  and  Vbc,  and  derives  the  Supply-side 
voltage  angle,  6s,  as  shown  in  Figure  3. 
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Figure  3.  Simulink  Model  of  Supply-side  Voltage  Angle  Calculation. 
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Since  the  SDC  only  has  four  voltage  inputs  it  is  important  to  be  effieient  about  the 
voltage  input  eonfiguration.  Further,  it  is  not  always  convenient  to  input  single  phase 
values  due  to  the  difficulty,  at  times,  to  tap  a  neutral  point.  This  is  why  it  is  suitable  to  use 
the  line-to-line  voltage  parameters  as  inputs.  The  derivation  of  the  Vab  and  Vbc  based 
voltage  angle  caleulation  as  described  in  [7]  and  [8]  is 


^ab  =  -^b’^bc  =^b-^c=^b  -{-^a-^bY 

Note  also  that  4+16+4=0  for  a  wye  conneeted  winding  set. 

Placing  into  matrix  form 
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'^bc 
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1 


-1 

2 


_  1 

2  1 

V/, 

~3 

-1  1 

We  define  the  matrix  in  dqO  terms  using 

cos(^)  +  Vj  eos^^-2^j  +  v^  eos^^  +  2^j 

(6»)  +  V,  cos  (<9  -  2^)  +  (-v^  -  V, )  eos  (<9  +  2^) 
s(^)-cos^^  +  2'^jj  +  v^  |cos^^-2'^j-cos^^  +  2'^jj 

and 

sin(^)  + Vj  sin^^-2^j  + sin^^  +  2^j 

sin(^)  +  Vj  sin  2.^1 +  (-v^  _y^^sin|^  + 2.^j 

^sin  {6)  -  sin  ^^  +  j  +  ^sin  j  -  sin  ^0  +  j 

After  combining  terms  the  line-to-line  voltage  equations  become 


(36) 


(37) 


(38) 


(39) 
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-2/ 

cos(6»J 

^ab 

“73 

sin(6>J 

-cos(e,  +  ^)_ 

ybc_ 

(40) 

Finally,  0s  is  calculated  using  the  ArcTan  funetion  with  Vq  and  Vd  as  the  inputs  for  the 
angle  caleulation.  The  ArcTan  function  gives  an  output  from  -Iti:  to  Itt  .  Sinee  Xilinx  is 
a  digital  bit-wise  based  program,  it  is  necessary  to  implement  the  last  block  of  Figure  3, 
which  converts  the  AreTan  output  to  a  wrapped  0  to  2^**  bit  output.  The  thetaconvl  bloek 
code  is  listed  in  the  Appendix. 

The  SDC  sends  the  sensed  power  supply  voltages  and  eurrents  and  processes  the 
data  in  order  to  send  the  eorrect  gate  signal  sequenee  to  the  SEMISTACK  BNC 
eonnections.  The  IGBTs  will  gate  on  and  off  and  send  the  supplied  voltage  to  an  internal 
eapaeitor  bank,  whieh  supplies  the  DC  Link  terminals.  Of  note,  eaeh  SEMISTACK  has 
it’s  own  eapaeitor  bank,  and  their  associated  external  terminals,  which  is  used  to 
physically  couple  the  Supply-side  and  DFIG-side  DC  link  to  form  a  eommon  Vdc-  The 
last  Supply-side  eireuit  input  is  Vdc  ,  which  is  sent  to  the  SDC  as  well.  This,  of  eourse,  is 
used  to  sense  and  maintain  DC  Link  voltage  at  200V.  It  is  also  used  in  a  protection 
feature  that  shuts  off  the  gate  signals  to  the  Supply-side  SEMISTACK  in  the  case  that  DC 
Eink  voltages  reaeh  240V.  The  implementation  of  this  proteetion  feature  utilizes  an  SR 
Flip-Flop. 

In  order  to  effectively  eontrol  the  DC  Link  voltage  an  inner  and  outer  PI 
controller  scheme  is  used.  The  outer  loop  controller  senses  Vdc,  as  compared  to  the  Vdc 
reference  of  200V,  and  generates  a  DC  reference  current  {iq_ref)  that  is  sent  to  the  inner 
loop.  The  inner  loop  is  controlled  in  the  qdO  synchronous  reference  frame  so  as  to  easily 
compare  iq  ref^o  the  measured  iq  in  DC  terms.  The  derived  6s  is  now  used  as  the  p  input, 
equation  (19),  for  the  qdO  equation  to  transform  the  supplied  currents  into  the 
synchronous  reference  frame.  Figure  4  shows  the  inner  and  outer  loop  Simulink  model. 


19 


Figure  4.  Inner  and  Outer  PI  Control  Loops  for  the  Supply-side  Controller. 

C.  DFIG-SIDE  CONTROLLER 

The  DFIG  rotor  is  powered  via  the  DFIG-side  SEMISTACK  IGBTs  through  the 
DC  Link.  The  DFIG  stator  is  powered  via  a  3-phase  AC,  1 20V,  60Hz  power  supply 
eorresponding  with  a  grid.  There  is  a  set  of  three  inductors  connected  in  series  to  the 
SEMISTACK,  downstream  of  the  DFIG-side  SDC,  which  decouples  the  SDC  sensors 
from  the  switching  action  of  the  IGBTs.  The  SDC  senses  the  stator  line-to-line  voltages, 
Vab  and  Vbc,  and  transforms  them  into  the  qdO  rotor  reference  frame  using  Oy  (from  the 
encoder)  as  [3  in  equation  (19).  With  the  stator  line-to-line  voltages  transformed  to  the 
dqO  rotor  reference  frame,  the  DFIG  slip  angle  and  frequency  can  be  derived  using  the 
same  method  as  described  in  the  Supply-side  controller  section  using  the  ArcTan  function 
with  Vq  and  Vd  as  the  inputs  for  the  angle  calculation.  The  DFIG  slip  angle  (dsu^  is 
important  since  it  will  be  used  to  control  the  DFIG  in  a  synchronously  rotating  dqO 
reference  frame,  with  the  g-axis  oriented  along  the  stator-flux  vector  position.  This,  in 
turn,  enables  decoupled  control  of  the  electrical  torque  {iqy)  and  rotor  excitation  current 
{Ui).  Traditionally,  OsUp  is  obtained  using  the  ArcTan  function  with  a  rotor  referenced  \(/qs 

and  \|/ds  as  inputs,  along  with  a  Yi  offset,  in  place  of  the  ArcTan  calculation  using  Vq  and 


Vd  as  the  inputs  for  the  angle  calculation.  The  offset  is  associated  with  the  relationship 
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between  current  and  flux  linkage  as  it  relates  to  voltage.  Since  the  OsHp  calculation  used  in 
this  thesis  is  based  on  voltage,  no  offset  is  needed.  The  voltage  based  calculation  of  OsHp 
requires  a  less  complex  algorithm,  and  thus  less  processing,  than  the  flux  linkage  based 
Osiip  calculation.  Both  experimental  and  simulated  tests  were  performed  to  compare  the 
two  Osiip  methods.  In  either  test  there  was  no  detectable  difference  in  the  results  for  this 
model’s  application.  Further  tests  should  be  performed  to  determine  the  extent  of 
constraint  between  the  two  methods  and  either’s  benefits.  Figure  5  shows  an  overview 
model  of  the  slip  angle  calculation.  The  internals  of  this  model  can  be  sufficiently 
described  by  reviewing  Figure  3. 


LINE-TO-LINE  qdO  TRANSFORMATION 
WITH  STATOR  VOLTAGES  IN  THE 
ROTOR  REFERENCE  FRAME 
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Figure  5.  Simulink  Model  of  DFIG-side  Slip  Angle  Calculation. 

In  addition  to  the  SDC  sensing  the  stator  line-to-line  voltages,  it  also  senses  the 
phase  A  and  B  rotor  currents  (phase  C  is  derived  using  phases  A  and  B).  The  rotor 
currents  are  utilized,  along  with  rotor  angular  frequency  (wr),  to  form  the  current  control 
loops  for  the  DFIG.  The  DFIG  PI  control  loops  consist  of  an  outer  loop  and  an  inner 
loop.  The  outer  loops  takes  Wr  as  an  input  and  compares  it  to  a  commanded  Wr_ref  ■  The 
resulting  difference  passes  through  the  speed  PI  controller  and  sends  a  reference  current 
{iqr  rej)  to  the  Current  control  loop,  or  inner  loop.  The  reference  current  passes  through  the 

current  PI  controller  and  compares  it  to  the  measured  iqr.  The  resulting  waveform  is 
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converted  to  a  per  unit  value  and  sent  to  the  SVM  block,  which  generates  the  gate  signal 
sequence  for  the  DFIG-side  SEMISTACK.  The  current  control  loop  is  used  to  command 
a  specified  Wr  as  controlled  via  iqr,  or  rotor  excitation  current  via  Ur-  Rotor  angular 
position  and  speed  are  derived  from  the  encoder,  which  will  be  discussed  later.  Figure  6 
shows  the  DFIG  controller  topology. 


Figure  6.  SIMULINK  DFIG  Controller  Topology. 

Both  the  Supply-side  and  DFIG-side  controllers  have  unique  differences  up  to  their 
respective  outputs,  however,  both  controller  outputs  go  to  the  exact  same  SVM 
architecture  in  order  to  control  their  commanded  parameters. 

D.  SPACE  VECTOR  MODULATION 

The  modulation  scheme  used  to  create  the  desired  3 -phase  waveforms  for  the 
Supply-side  and  DFIG-side  circuits  is  SVM.  SVM  is  a  specific  form  of  Pulse  Width 
Modulation  (PWM)  in  which  an  algorithm  involving  space  vectors  are  used  to  control  the 
on  and  off  times  of  pulsed  signals.  The  generated  signals  then  drive  the  IGBT  gate  signals 
to  the  SEMISTACK  allowing  the  user  to  create  a  waveform  of  any  magnitude  and 
frequency  desired.  The  SVM  approach  utilized  in  this  controller  is  obtained  from  [8-9]. 

The  SEMISTACK  produces  a  3 -phase  output  voltage,  as  supplied  from  the  DC 
Link,  that  is  controlled  by  turning  on  and  off  its  six  IGBTs.  Figure  7  shows  both  the 
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SEMISTACK-IGBT  configuration  and  SVM  Hexagon  for  which  the  on  and  off 
switching  states  are  based  on.  The  gJ-axis  is  overlaid  on  the  SVM  Hexagon  and 

represents  the  axis  of  the  stationary  reference  frame  variables  y^andv^^,  which  are 
shown  entering  the  SVM  block.  The  magnitude  of  the  vectors  and  is  equal  to  the 
vector  V*  and  6  is  equal  to  the  ArcTan  of  the  two  vectors. 


p  bus 


V. 
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f 

V, 

1 

3^  ri 

^ . i 

V.  - 

Rotor 


Figure  7.  SEMISTACK-IGBTs  and  SVM  Hexagon  [From  8]. 


There  are  six  sectors  on  the  hexagon  and  eight  possible  on  and  off  states  that  are  used  to 
produce  the  V*  vector.  The  vectors  V\  and  V2,  when  summed,  are  also  equal  to  the  vector 
V*.  In  Sector  I,  V\  and  V2  correspond  to  the  states  (/>,«,«)  and  (/>,/),«) .  The  /?and  n 

correspond  to  the  positive  or  negative  IGBT  bus  signals  for  the  respective  phases V^,V^, 
andF^.  The  zero  vectors,  which  are  represented  by  the  states  (p,p,p)  and  («,«,«),  are 

not  shown  but  are  represented  as  vectors  into  and  out  of  the  page.  The  magnitude  of  the 
vectors  V\  and  V2  correspond  to  the  amount  of  time  spent  on  the  switching  states  in  each 
sector.  The  duty  cycles  T^  and  are  the  times  spent  on  each  cycle  for  the  vectors  Fj 


and  .  The  total  time  spent  for  one  switching  period  is 


K  = 


2TV 

d 


dc 


3T 


(41) 
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(42) 


2TV 


3T 


The  law  of  sines  on  any  one  of  the  seetors  in  SVM  Hexagon  will  produee 

2F*  _  V,  _ 

sin  (6O‘’-0)  sin(^) 

(43) 

and  substituting  equations  (41)  and  (42)  into  (43)  is  used  to  find  the  time  of  the  duty 
eycles  for  the  vectors  and  F2 


2;=-^z;sin(60"-«) 

^  dc 

T,  =  '^T,sm(e) 


T  =T  +T +T 


(44) 

(45) 

(46) 


SVM  has  the  ability  to  minimize  switching  loss  and  harmonic  distortion  by  controlling 
the  order  in  which  the  states  are  applied  so  that  only  one  gate  is  fired  between  sectors  [8- 
9],  The  switching  pattern  for  each  state  and  sector  used  in  this  thesis  is  shown  in  Table  1. 
The  SVM  algorithm  is  implemented  inside  each  SDC  where  it  is  oversampled  to 
converge  on  the  performance  of  a  true  continuous  sinusoidal  signal.  The  digital 
implementation  of  the  SVM  scheme  used  for  this  thesis  is  discussed  in  detail  in  [8-9]  and 
the  Simulink  diagrams  can  be  found  in  Appendix  C.  Figure  8  is  a  block  diagram  of  the 
functions  for  each  process  in  the  SVM  block. 
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Table  1 .  Space  Vector  Modulation  Switching  Pattern  [From  9], 


Figure  8.  SVM  Digital  Implementation  Diagram  [From  8]. 
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E.  ENCODER  IMPLEMENTATION 


The  encoder  used  for  the  DFIG  rotor  position  and  speed  input  is  the  Microtech 
Laboratory  Inc.  (MTL)  MES-20-200P  C.  Its  implementation  is  discussed  in  [10], 
however  several  modifications  were  made  to  commission  the  encoder  in  this  thesis.  The 
encoder  uses  a  single  shaft,  200  pulse,  open  collector  output  design  to  produce  an  Ak,  Bk 
and  Zk  pulse  sequence  as  seen  in  Appendix  A  and  Figure  9  (experimental  Chipscope 
capture). 


Figure  9.  Experimental:  Encoder  A,  B  and  Z  Pulse  Configuration. 

Of  note,  the  imperfections  shown  in  Figure  9,  and  some  subsequent  figures,  are 
caused  by  a  Chipscope  plot  function  artifact  due  to  Chipscope  timing  discrepancies. 
Pulses  Ak  and  Bk  are  in  quadrature,  and  the  Zk  pulse  occurs  once  every  revolution  (200 
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pulses)  in  alignment  with  the  B  pulse  period.  In  order  to  derive  the  rotor  position  from 
this  encoder  input  a  Simulink  block  was  created  to  use  the  Ak  and  Bk  pulses  as  an 
incremental  counter  and  the  Zk  pulse  as  a  one  revolution  reset.  Figure  10  is  the  Simulink 
block  implemented  for  this  function. 


Figure  10.  Architecture  Used  to  Derive  Rotor  Angular  Position  from  the  Encoder. 

The  concatenation  block  takes  Ak,  Ak-i  and  Bk  and  forms  a  three  bit  word  and 
relates  it  to  a  two  bit  word  (Pk  and  Nk  combined).  The  two  bit  word  is  compared  to  the 
predetermined  encoder  vector  in  the  ROM  block  (Pk  and  Nk  is  also  used  to  enable  the 
counter).  The  predetermined  encoder  vector  is  specific  to  the  encoder  characteristics  and 
is  used  to  indicate  whether  the  encoder  shaft  is  turning  clockwise  (binary  1  0)  or  counter¬ 
clockwise  (binary  1  1),  hence  count  up  or  down.  A  rising  and  falling  edge  architecture  is 
chosen  for  the  encoder  algorithm,  which  effectively  doubles  the  encoder  resolution  from 
200  triggers  per  revolution  to  400  triggers  per  revolution.  The  two  bit  word  is  then  passed 
through  the  Simulink  Slice  block,  which  separates  the  word  into  two  one  bit  words.  Pk 
and  Nk  trigger  the  counter  to  count  up  or  down  with  a  range  of  0  to  2'*^  bits  with  a  wrap 
feature  triggered  by  the  reset  pusle.  Table  2  shows  a  detailed  account  of  the  encoder 
algorithm. 
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Ak 

Ak-1 

Bk 

Pk  Nk 

(Logical  OR) 

Enable  Up 

Counter  State 

Encoder  Vector 

0 

0 

0 

N/A 

N/A 

N/A 

N/A 

0 

0 

1 

N/A 

N/A 

N/A 

N/A 

0 

1 

0 

1  0 

1  1 

Count  down 

Falling  edge  CCW  (2) 

0 

1 

1 

0  1 

1  0 

Count  up 

Falling  edge  CW  (1) 

1 

0 

0 

0  1 

1  0 

Count  up 

Rising  edge  CW  (1) 

1 

0 

1 

1  0 

1  1 

Count  down 

Rising  edge  CCW  (2) 

1 

1 

0 

N/A 

N/A 

N/A 

N/A 

1 

1 

1 

N/A 

N/A 

N/A 

N/A 

Table  2.  Encoder  Truth  Table  with  Directed  Actions. 


The  Zk  pulse  provides  inputs  to  an  AND  gate  as  seen  in  Figure  10.  Since  it  occurs  only 
once  per  revolution,  it  is  straightforward  to  see  why  it  is  used  for  the  counter  reset.  The 
counter  output  is  a  true  account  of  the  encoder  position,  however  the  actual  rotor 
electrical  position  is  not  accurate  until  the  encoder  zero  point  is  calibrated  to  be  aligned 
with  the  rotor  electrical  zero  point. 

The  encoder  calibration  procedure  began  by  a  detailed  study  of  the  rotors  physical 
winding  layout.  The  DFIG  is  a  4-pole,  3-phase  induction  machine  with  24  slots  and  no 
accessible  rotor  windings  neutral  point.  A  simple  method  of  aligning  the  rotor  to  a  zero 
point  electrical  position  would  be  to  supply  power  to  the  phase  A  windings  of  the  rotor 
and  stator  and  allow  the  magnetizing  flux  to  align  the  two.  Without  access  to  the  rotor 
winding’s  neutral  point,  this  method  became  moot.  The  second,  and  possibly  more 
accurate,  method  to  align  the  rotor  to  its  zero  point  electrical  position  was  to  supply 
power  to  the  rotor  and  stator  phases  via  line-to-line.  After  several  experiments  with  the 
line-to-line  method,  there  were  multiple  convergences  to  a  rotor  zero  point  electrical 
position  to  be  3 1 8  triggers  out  of  400.  This  offset  was  added  to  the  encoder  algorithm  as 
the  calibration  point.  Finally,  the  calibrated  encoder  signal  was  passed  through  a  series  of 
conversions  so  that  it  corresponded  to  the  rotor  electrical  angular  position  and  also  had  a 
wrapping  range  of  0-2'*^  bits  per  revolution.  Figure  11  shows  the  experimental  results  of 
the  calibrated  rotor  electrical  angular  position. 
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Figure  1 1 .  Experimental:  Calibrated  Rotor  Electrieal  Angular  Position. 

The  last  input  the  encoder  provided  was  the  electrical  angular  frequency,  Wr.  This 
derivation  was  simply  a  differentiation  of  rotor  electrical  angular  position  with  respect  to 
time.  Due  to  the  limitations  of  Simulink’s  iterative  features,  the  output  of  this 
differentiation  was  crude  and  noisy.  However,  after  passing  this  output  through  a  first 
order  low  pass  filter  with  a  comer  frequency  of  lOHz,  the  signal  was  surprisingly 
accurate  and  proved  to  be  very  functional  for  this  application. 


29 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


30 


IV.  RESULTS 


A.  OVERVIEW 

In  this  thesis,  we  emulate  a  Wind  Turbine  Generator  by  driving  a  DFIG  via  a  DC 
motor  with  variable  input  torque  capability.  As  discussed  in  [1-4],  we  can  use  a  DFIG 
with  back-to-back  SVMs  to  accomplish  a  decoupled  Supply-side,  and  DFIG-side,  control 
scheme  while  allowing  power  flow  to  occur  between  the  two  systems. 

The  two  circuits  of  concern  are  the  DFIG  and  Supply-side  circuits.  The  DFIG  and 
Supply-side  circuits  are  electrically  coupled  with  back-to-back  (DFIG-side  and  Supply- 
side)  SVMs  that  are  coupled  through  a  DC  Link  consisting  of  a  capacitor  bank.  The  back- 
to-back  SVMs,  with  DC  Link,  provide  bi-directional  power  flow  between  the  DFIG  rotor 
and  Supply-side  power  supply.  Bi-directional  power  flow  is  achieved  when  the  DFIG  is 
controlled  at  supersynchronous  speeds,  and  input  torque  is  enough  to  overcome  losses, 
and  subsynchronous  speeds  (Supply-side  provides  power  to  rotor). 

B.  SUPPLY-SIDE  EXPERIMENTS 

The  Supply-side  power  supply  (3 -phase  AC)  provides  power  to  the  Supply-side 
circuit  (stepped  down  by  a  VARIAC)  and  the  Stator  of  the  DFIG.  The  Supply-side  circuit 
senses  the  supply  side  voltage,  current  and  DC  Link  voltage  and  sends  the  sensed  inputs 
to  a  FPGA.  The  FPGA  is  programmed  via  a  Simulation  software  package  called  Simulink 
to  achieve  the  desired  control  of  the  specified  parameters.  For  the  Supply-side  circuit  we 
use  the  mentioned  sensed  inputs  to  derive  control  of  Vdc  so  that  it  is  maintained  at  200V 
via  SVM.  There  is  also  a  built-in  protection  feature  in  the  Supply-side  circuit  that  will 
turn  off  the  SVM  if  a  VDC  of  240V  is  reached.  Within  the  Supply-side  controller 
architecture,  however,  there  are  parameters  that  can  be  used  to  control  the  system  power 
factor.  The  research  of  this  feature,  and  its  effects  on  the  grid,  was  beyond  the  scope  of 
this  thesis,  but  some  simulation  and  experimental  tests  were  conducted  to  show  the 
performance  of  the  controller.  If  a  system,  connected  to  the  grid,  has  explicit  control  of  its 
power  factor,  it  can  be  used  to  reduce  transmission  line  losses  by  controlling  the  amount 
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of  reactive  power  being  distributed.  Figures  12  and  13  show  the  simulation  and 
experimental  steady  state  shifts  in  power  factor  when  id_ref  was  commanded  at  2A  and  - 
2A. 


Figure  12.  Simulation:  Supply-side  Voltage  and  Current  with  id_refSit  2 A  and  -2A. 
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Figure  13.  Experimental:  Supply-side  Voltage  and  Current  with  id  ref^^  2A  and  -2A. 


C.  DFIG-SIDE  EXPERIMENTS 

The  DFIG  circuit  power  supply  is  provided  by  the  DC  Link  and  DFIG-side  SVM. 

The  DFIG  circuit  senses  rotor  current,  stator  voltage,  rotor  speed,  and  rotor  electrical 

position  (as  provided  by  an  encoder).  These  inputs,  like  the  Supply-side  circuit,  are  then 
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sent  to  an  FPGA,  which  is  programmed  to  control  the  DFIG  parameters.  The  DFIG 
circuit  will  take  the  rotor  speed  input  and  derive  a  current  reference  that  is  used  to  control 
the  rotor  speed  using  a  PI  controller.  It  is  important  to  show  that  we  can  control  the  DFIG 
rotor  speed  since  traditional  Wind  Turbine  Generators  would  use  turbine  blade  pitch 
control  in  tandem  with  a  synchronous  generator  to  provide  a  constant  speed,  and  hence  a 
constant  output  frequency.  Further,  it  is  noteworthy  that  a  wind-speed  profde  can  be 
adapted  to  a  specific  Wind  Turbine  (outfitted  with  a  DFIG)  in  order  to  program  the  DFIG 
to  run  at  specified  speeds  for  a  given  wind  speed.  This  tool  can  be  used  to  provide  an 
optimized  rotor  speed  for  a  given  wind  speed.  For  our  application  we  will  show  that  the 
DFIG  rotor  speed  can  be  controlled  (subsynchronous,  synchronous  or  supersynchronous) 
for  any  given  input  torque  and  thus  distributes  the  power  generated  by  the  DFIG  as 
predicted  between  the  rotor  and  stator.  Further,  this  experiment  will  demonstrate  the  bi¬ 
directional  power  flow  of  the  back-to-back  SVMs  showing  power  recovery  via  the  DFIG 
rotor. 


Figure  14.  Simplified  Circuit  of  Wind  Turbine  DFIG  System. 
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As  illustrated  in  Figure  14,  the  Supply-side  SVM  converts  a  3-phase  AC  source  to 
a  constant  DC  Link  voltage  of  200V.  The  DC  Link  voltage  is  then  inverted  into  a  3 -phase 
waveform  via  the  DFIG-Side  SVM  (during  modes  for  which  power  is  being  delivered  to 
the  DFIG).  The  digital  control  for  the  Insulated  Gate  Bipolar  Transistor  (IGBT)  network 
is  implemented  via  the  FPGA.  The  FPGA  interface  samples  the  rotor  current  and  the 
stator  voltage  of  the  DFIG  through  an  A/D  converter.  There  is  an  input  from  the  encoder, 
Or,  to  provide  the  position  and  speed  of  the  rotor  to  the  FPGA  through  the  interface  circuit 
board.  An  algorithm  then  calculates  the  slip  frequency  of  the  DFIG,  which  is  used  to 
derive  the  needed  rotor  current  frequency  and  amplitude  to  maintain  a  specified  rotor 
speed,  and  therefore  a  controllable  output  frequency.  Since  the  Supply-side  and  DFIG- 
side  controller  topologies  are  nearly  identical,  it  follows  that  the  DFIG-side  controller  can 
control  system  power  factor  in  the  same  way  that  the  Supply-side  did.  Figures  15  and  16 
show  the  simulation  and  experimental  steady  state  shifts  in  power  factor  when  id  ref  was 
commanded  at  2A  and  -2A. 
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Figure  15.  Simulation:  DFIG-side  Voltage  and  Current  with  id_ref^i  2A  and  -2 A. 
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Figure  16.  Experimental:  DFIG-side  Voltage  and  Current  with  id_ref^i  2A  and  -2 A. 

Figure  17  demonstrates  an  input  torque  transient  representing  a  large  wind  gust. 
By  commanding  a  specified  supersynchronous  speed,  while  taking  into  account  windage 
and  friction  losses  of  the  system,  it  is  clear  that  the  input  torque  transient  produces  a 
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power  flow  reversal  in  the  rotor.  In  other  words,  the  rotor  mode  transitions  from  power 
flow  from  the  rotor  (negative  rotor  current)  to  power  flow  to  the  rotor  (positive  rotor 
current). 


3 
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Figure  17.  Simulation:  Input  Torque  Transient  Showing  Rotor  Power  Reversal. 


Figures  18  and  19  (experimental  results)  parallels  this  transition  with  the  steady 
state  captures  corresponding  to  the  simulation  transient.  The  experimental  results  were 
run  through  a  low  pass  filter  with  a  corner  frequency  of  80  Hz  in  order  to  create  a  useful 
plot  of  the  event.  Due  to  windage  and  friction  losses,  a  commanded  supersynchronous 
speed  was  ordered  to  achieve  the  power  flow  transition  results  during  the  constant  rotor 
speed  experiment. 
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Figure  18.  Experimental:  Rotor  Current  and  Voltage  with  Power  Flow  from  the  Rotor. 
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Additionally,  Figure  17  (simulation  model)  shows  a  smooth  transition  as  the 
power  flow  reverses  between  power  flow  from  the  rotor  to  power  flow  to  the  rotor.  This 
transition  is  modeled  using  a  constant  commanded  rotor  speed  along  with  a  10%  decrease 
in  input  torque. 

Figure  20  demonstrates  a  large,  constant  input  torque  representing  a  constant  wind 
speed.  By  increasing  the  commanded  rotor  speed  by  30%,  while  taking  into  account 
windage  and  friction  losses  of  the  system,  it  is  clear  that  the  rotor  speed  transient 
produces  a  power  flow  reversal  in  the  rotor  that  corresponds  to  subsynchronous, 
synchronous  and  supersynchronous  speeds.  In  other  words,  the  rotor  mode  transitions 
from  power  flow  to  the  rotor  (positive  rotor  current)  to  power  flow  from  the  rotor 
(negative  rotor  current). 
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Figure  20.  Simulation:  Increase  in  Rotor  Speed  with  a  Constant  Input  Torque. 

Figures  21  through  23  (experimental  results)  parallels  Figure  20  with  the  steady 
state  captures  corresponding  to  the  simulation  transient.  The  experimental  results  were 
run  through  a  low  pass  filter  with  a  corner  frequency  of  80  Hz  in  order  to  create  a  useful 
plot  of  the  event. 
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Figure  2 1 .  Experimental:  Rotor  Current  and  Voltage  While  Subsynchronous. 
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Figure  22.  Experimental:  Rotor  Current  and  Voltage  While  Synchronous. 
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Figure  23.  Experimental:  Rotor  Current  and  Voltage  While  Supersynehronous. 


In  eonclusion,  we  showed  that  we  ean  control  the  DFIG  rotor  excitation  to 
provide  a  constant  output  frequency.  Our  application  demonstrated  that  the  DFIG  rotor 
speed  can  be  controlled  (subsynchronous,  synchronous  or  supersynehronous)  for  a  given 
input  torque  and  thus  distributes  the  power  generated  by  the  DFIG  as  predicted  between 
the  rotor  and  stator.  Finally,  this  experiment  also  demonstrates  the  bi-directional  power 
flow  of  the  back-to-back  S  VMs  showing  power  recovery  via  the  DFIG  rotor. 
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V.  CONCLUSIONS  AND  FUTURE  RESEARCH 


A.  CONCLUSIONS 

The  Wind  Turbine  Power  Generation  Emulation  via  Doubly  Fed  Induction 
Generator  Control  System  was  designed,  simulated,  constructed  and  tested  using  standard 
engineering  principles  and  theory.  Since  the  research  and  simulation  phases  were  aptly 
utilized,  the  development  and  construction  phases  were  relatively  inexpensive  and  short. 
A  series  of  test  were  conducted  and  validated  the  design.  The  experimental  results 
showed  that  the  system  responded  according  to  the  simulation. 

With  a  baseline  prototype  that  functions  correctly,  the  Naval  Postgraduate  School 
now  has  a  means  to  build  a  large  scale  prototype  for  laboratory  purposes  and  future 
research.  The  future  research  on  this  topic  will  be  above  and  beyond  current  publications, 
and  will  have  an  enabled  freedom  to  exploit  the  abilities  of  the  DFIG. 

B.  FUTURE  RESEARCH 

A  third  SDC  should  be  used  to  take  control  of  the  DC  motor.  This  Wind-side  SDC 
can  be  programmed  to  mirror  current  wind  speed  profdes  that  are  available  from  the 
National  Renewable  Energy  Laboratory.  With  a  wind  speed  profde  controlling  the  input 
torque  of  the  DFIG,  the  DFIG-side  SDC  can  be  slightly  modified  to  include  a  rotor  speed 
lookup  table.  This  lookup  table  can  be  used  to  derive  an  optimal  rotor  speed  for  a  given 
wind  speed,  and  thus  increase  the  efficiency  of  DFIG  power  extraction. 

Second,  DFIG-side  feed  forward  terms,  based  on  slip  speed  and  flux  linkages,  can 
be  implemented  as  described  in  [1].  The  feed  forward  terms  should  lower  the  DFIG-side 
PI  controller  time  constant  without  any  detrimental  effects  on  stability. 

Finally,  Supply-side  feed  forward  terms,  based  on  the  stators  synchronous  angular 
frequency  and  stator  winding  inductance,  can  be  added.  Again,  these  terms  should  lower 
Supply-side  PI  controller  time  constants  without  any  detrimental  effects  on  stability. 
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APPENDIX  A:  DATASHEETS 


SEMISTACK  -  IGBT 


Three-phase  rectifier  + 
inverter  with  brake 
chopper 

SEMITEACH  -  IGBT 
SKM  50  GB123D 
SKD  51 

P3/250F 

Features 


Circuit 

Inre  (A) 

Vac  /  Vcjmiajj 

Types 

BSCI 

30 

440  y  750 

SEMiTEACH  -  IGBT 

Symbol 

Conditions 

Values 

Units 

'Irma 

nooveiioari 

at} 

A 

tGBT  -  4x  SKM  50  GB  1230 

VcES 

1200 

V 

VcEtSAp 

VaES 

k=  SOA.  VGf=  15V.  chip  level;  Tj=  25{12S)^C 

2.7  (3.5} 
±20 

V 

V 

(c 

ktt 

Tcib^  25 1  sore 

T™=25  mrC:  tp=1ms 

50(40} 

100(80} 

A 

A 

VhJmaiJ 

HeeSfrer  -  1k  SKD  5iyi4 

without  filter 

3*430 

V 

with  filter 

3x330 

V 

VdCttik 

DC  Capaotor  bank  -  Electrolytic  2x  2200tJFi'4[X}V 

total  erjulvalenl  capacitance 

max.  DC  voltage  applied  to  Ihe  capacilor  bank 

1100/800 

750 

pF/V 

V 

Dover  -  4x  SKHJ  22 

Power 

supply 

Current 

oy  15 

V 

consump 

lion 

max;  per  driver 

16 

mA 

TtiermaJ 

irip 

Normally  Open  lype  (NO) 

71 

•  Muhti-function  IGBT  converter 

'  Transparent  enclosure  to  allow 
visualization  of  every  part 

•  IP2k  protection  to  minirriize 
safety  hazards 

'  Ejftemal  banana/BNC  type 
conr>ectors  for  all  devices 

•  Integrated  drive  unit  offering 
short-circuit  deleciiorv'cut-off. 
power  supply  failure  detection, 
^teriock  of  IGBTTs  -i-  galvanic 
isolatiosn  of  the  user 

•  Forced-air  cooled  heatsink 

Typical  Applications 

•  Education:  One  stack  can 
simulate  almost  all  existing 
industrial  applications: 

-  3-phase  inverter+brake  chopper 

-  Buck  or  boost  converter 

-  Single  phase  inverter 

-  Single  or  S-phase  rectifier 

^1  Pholcj  nan-CDfrtrartual 


n  ri  _  n  _n  _n  n  n 


Ql 

nanenBris 

rtf 
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'VfV: 
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IC 

■  Mfl-  1 

P  ow 
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o 

r  -iUj  ■ 

• - 3oa - ■ 

General  dimensions 


TNs  tGctifiical  Information  specrfica  aflmacorxtuciDr  devices  but  promises  no 
charact-gristics.  Mo  warranty  or  guarantee  expressed  or  implied  is  made  regarding 
delivery,  perfcfmance  or  suitabiFity. 
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Model  8231  Three-Phase  Wound-Rotor  Induction 
Motor 


Each  phase  of  the  stator  windings  of  this  motor  is  inde¬ 
pendently  terminated  and  identified  on  the  faceplate  to 
permit  operation  in  either  deita  or  star  (wye)  configura¬ 
tion.  The  rotor  windings  are  brought  out  to  the  faceplate 
viaerfemal  slip  rings  and  brushes.  This  nrachine  can  l>e 
used  as  a  wound-rotor  induction  motor,  phase  shifter, 
single-phase  variable  coupling  transtofmer,  three-phase 
transformer,  selsyn  control,  tequency  converter  or 
asynchronous  induction  genefator.  The  speed  of  this 
machine  can  i>e  controlled  through  the  use  of  the  Three- 
Phase  Rheostat  (Model  8731). 

SPECIFICATIONS 


Model  8211  DC  Motof/Generahw 

l2<W20BV-60  Hz 

220raBOV-  50  Hz 

240M15V-50HZ 

Power  Requirwnent 

12a,'2DBV 

220.^  V 

24[1.'^15V 

Rating  Motor  Output  Power 

Oeft«ator  Output  PowEf 

Armature  Voltage 

Shud  Held  Voltage 

PuR  Load  S^eed 

Ful  Load  Motor  Current 

Pull  Load  C5eneratof  Current 

175  W 

120  W 

now 

12DW 

120V-DC 

230  V- DC 

240  V- DC 

120V-DC 

220V~  DC 

240  V- DC 

1GDD  r^in 

1503  nhm 

1 500  r/mfi 

ZB  A 

1JA 

1.1  A 

1  A 

D.5A 

D.5A 

Physical  Char^oteristios  Oim^nrsions  (N  sc  W  k  D) 

Net  Weight 

30flK2B1  x+40rT¥n{12.1  x  11.5 x  17.3  h) 

14.1  kg  (31  fc) 

Model  8221  Four-Pole  SquiiTel-Cag«  [nduction  Motor 

12[M208Vwea  Hz 

22ftraB0V^  50  Hz 

240/415  V~.  50  Hz 

Power  Requrement 

12D™V 

220i'3S0V 

24D.'415V 

Ratng  Output  Powe^ 

Stator  Voltage 

Pull  Load  S^peed 

Pull  Load  Cuirent 

175  vy 

i2a-2aa  v. 

22ai'3BD  V,  3-phase 

240^415  V.  3-phase 

167D  nhiin 

13flOD'mii 

1 395  r/mr 

1ZA 

0.52  A 

0.46  A 

Physcal  Characteristics  DintefMions  (H  jc  W  x  □) 

Net  Weight 

305x261  X  440  rrmf  12..1  x  11.5x  17J  in> 

13.S  kg  (29  J  lb| 

Model  8231  Three-Phase  Wound-RolDr  Induction  Motor 

I2ft!208  V-  M  Hz 

221W3B0  V-  50  Hz 

240/415  V-  50  Hz 

Power  Requrement 

12t].'2D5  V 

220.^  V 

24Q/415V 

Ratrg  Output  Fnwef 

Stior  Voltage 

Rotor  Voltage 

Pull  Load 

Full  Load  Current 

175W 

12fl'2aa  V.  3rphasE 

220i'3BD  V,  a^jhase 

24E1.415V.  3-phase 

efliT04V.  3-phase 

llQi'IBDV.  3-phase 

12[L'3a0V.  3-phase 

1 50D  r^in 

1240  rfrrun 

1315  r.'rrwi 

1.3  A 

0.53  A 

D.4BA 

Physical  CharadensScs  DintEniSians  (N  x  W  x  D) 

Net  Weight 

305  k291  X  440  rrwi  02-1  k  11.5k  17.3  in) 

14  kg  (BD.a  fc) 
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POWER  SUPPLY 
MODEL  8821 


The  Power  Supply  provides  fixed  and  variabte  AC 
and  DC  voltage  sources,  all  terminated  by  color-coded 
4  mm  safety  sockets.  Independent  circuit  breakers,  reset 
at  the  front  panel,  protect  the  input  to  and  output  from  the 
Power  Supply.  Indicator  lamps  monitor  tfre  presence  of 
input  voltage  in  each  phase.  When  a  phase  leg  of  the 
site's  power  service  is  out,  the  lamp  goes  off  to  reflect  this 
condition. 


A  voltmeter,  connected  through  a  selector  switch, 
monitors  the  variable  AC  and  DC  outputs  and  fixed  DC 
output  A  24  V  AC  output  provides  a  low-voltage  supply 
required  to  operate  othe  r  EI^S  equ  ipm  ent  such  as  meter¬ 
ing  modules  and  modules  used  in  the  Power  Electronics 
Training  System. 


SPECIFiCATIONS 


Mo[k-l  0021  -  Powgr  Suppriy 

120/208  V- 60  Hz 

220/380  V  -  50  Hz 

24*^415  V- 50  Hz 

Inpvt  UneVDltagg 

Line  Current 

Service  IrtEtallatiDn 

12a.'2D6V 

22a.'30a  V 

24a'415V 

16  A 

10  A 

20  A.  3-phase  .  5  wires,  star  (WyeJ-  oonneded.  hduding  neutral  and 
gro^TC 

Outputs  Faed  AC  S-Phase 

Variable  AC  3-Phase 

Variable  DC 

Faed  DC 

Low  Power  AC 

i2U.{20^V-  16  A 

220/350  V- 10  A 

24D.'415V  -  IDA 

0-120.^00  V-5A 

D-22a™V-3A 

0-240/416  V-3  A 

D-12D  V-OA 

D-220  V  -  5  A 

0-240  V-6  A 

12D  V-2A 

220  V  -  1  A 

240V-1  A 

24  V-3A 

WaS  Outlet  induded 

NEMAL21-20 

NEMA  L22-2D 

CLIPSAL  56SO520 

Power  Cord 

3  m  ( 1 D  ft  j 

Physical  Charact&^titjs  Dimensofts  (H  k  W  x  □) 

Net  WeigH 

303  )t  237  K  5DD  rwn  (12.1  )c1 1.3  x  10.7  in) 

13.4  kg  (40.5  lb) 
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APPENDIX  B:  MATLAB  M-FILES 


A.  MATLAB  INITIAL  CONDITIONS  FILE 
DFIG.m 

Vdc=200;  %%%DC  Link  Voltage  Setpoint 

V_phase=sqrt (2 ) ^120 ;  %%%Supply-side  voltage 

Ks=2/3^[1  -1/2  -l/2;0  -sqrt(3)/2  sqrt ( 3 ) /2 ; 1/2  1/2  1/2];  %abc  to  qdO 

transformation  in  the  stationary  frame 

f_fund  =  60;  %%%base  frequency 

omega_b  =  2'^pi'^60;  %%%wr 

oversample=2 ;  %%%VSM  oversample 

step_ct=l ; 

pulsect=1800/step_ct ; 
clock_f  req=25e6  ; 
tstep  =  40e-9^step_ct ; 
tstop=4 ; 

o.o.o.o.o.o.o.o.o.o.o.o.o.o.o.o.o.C-hrD-l-rN-r'  DT  r'rD-i  ino9^9^9^9^9^9^9^9^9^9^9^9^9^9^9^9^9^9^9^9^9^9^ 
oooooooooooooooooOLdLU-L  irJ-  kirci-Lliooooooooooooooooooooooo 

Kp_vdc= . 1 ; 

Ki_vdc=10 ; 

Kp_i_s=10  ; 

Ki~i~s=20; 

Ki_i_s_icq=0 ; 

Ki_i_s_icd=0 ; 

9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9' 
oooooooooooooooooooooooooooooooooooooooooooooooooooooo 


%%%clock  frequency 
%%%step  size 


Lf=300e-6;  %%%choke  inductor 

Rf=.05; 


Amat_indl 
Bmat_indl 
Cmat_indl 
Dmat  indl 


-inv([Lf  -Lf;Lf  Lf +Lf ] ) ^Rf ^ [ 1  -1;1  2]; 
-inv([Lf  -Lf;Lf  Lf+Lf ] ) ; 

[1  0  ;0  1  ;-l  -1  ];  %Ic  =  -la-Ib 
zeros (3,2) ; 


one_zero_state=0 ;  %Set  to  one  so  that  only  one  zero  state  is  used 

in  modulation 
if  one_zero_state  ==  1 
gainl  =  1; 
gain2  =  0; 

else 

gainl  =  1/2; 
gain2  =  1; 

end 


twopiby3  =  2'^pi/3; 
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poles  =  4; 

polesby2 J=poles/2/ ( . 089) ; 

base  torque  calculation 

Tno9'9'9'9'9'9'9'9'9'9-9'9'9' 
ocL-Lllo  ooooooooooooo 

Kp_i=80; 

Ki  i=100; 


TL=200/ (2/poles*omega_b) ;  %%%DFIG 
%%%%%%%%%%%%%Rotor  Speed/Current  PI 


Kpgain_speed= . 0333 ; 

Kigain_speed= . 0033; 

oooooooooooooooooooooooooooooooooooooooooooooooooooooo 


rs=12 ; 
rr  =  4  ; 
Xls  =9; 
Xm  =180; 
Xlr  =  9; 


%%%stator  res 
%%%rotor  res 


D= (Xls+Xm) ^ (Xlr+Xm) -Xm^2 ; 
rsbyXls  =  rs/Xls; 
rrbyXlr  =  rr/Xlr; 

Xaq  =  1/ (1/Xm+l/Xls+l/Xlr) ; 

Xad  =  Xaq; 

XaqbyXls  =  Xaq/Xls; 

XaqbyXlr  =  Xaq/Xlr; 

XadbyXls  =  Xad/Xls; 

XadbyXlr  =  Xad/Xlr; 

XaqbyXm  =  Xaq/Xm; 

XadbyXm  =  Xad/Xm; 
psi_qsic=0 ; 
psi_dsic=0 ; 
psi_qric=0 ; 
psi_dric=0 ; 
omegar_ic  =  omega_b; 

Kp_iqs=50 ; 

Ki_iqs= . 1 ; 

Kp_iqr=50 ; 

Ki_iqr= . 1 ; 

F_mat  =[0001; 1120; 223 
0  mat  =  F  mat; 


0 ;  3  3  0  0  ]  ; 


ooooooooooooooo 


o  o  o 
o  o  o 


Encoder 


Speed 


\T(=^  n-\-r^-r9-9'9-9'9'9'9-9'9-9'9-9'9-9-9-9-9-9'9-9'9'9'9-9'9-9'9-9'9'9'9-9'9-9'9-9'9- 
VtdOLUJLooooooooooooooooooooooooooooooooooooo 


lndex=[l:2^12]  ; 
reciprocal=2^-6 . / Index; 


2'9'2'9'2'2'2'9'9'2'2'9'2'9'2' 

ooooooooooooooo 


o  o  o 
o  o  o 


Encoder 


Vector 


2'2'2'9'2'2'2'9' 

oooooooo 


o, 

o 


output_vec= [ 0 ; . . . 


0 

2 

1 


%Motor  is  turning  in 
%Motor  is  turning  in 


the  CCW  falling  edge 
the  CW  falling  edge 
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1; . . . %Motor 
2; . . . %Motor 

0 ;  .  .  . 

0]  ; 

ooooooooooooooooooooooo 

o  o  o  o 
o  o  o  o 


is 

is 


turning  in  the 
turning  in  the 


CW  rising  edge 
CCW  rising  edge 


oooooooooooooooooo 


ooooooooooooooooooooooooooooo 
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B. 


MATLAB  M-FILE  USED  FOR  SPACE  VECTOR  MODUALTION 


cverflow3.m 

function  [sectorl,  sector2,  sectorS,  sector4,  sectors, 
overflows (x) 

%gain  =  xf ix ( {xlUnsigned, 10, 7 } , 2 . 359296/3) ; %for  60  hz 
gain  =  xf ix ( {xlUnsigned, 10, 7 } , 2 . 359296) ; %for  180  hz 
%  t  emp  v= g  a  i  n  X  ; 
tempv=x; 
if  tempv<=17 1-1 

sectorl=xfix ( {xlBoolean} , 1) ; 
sector2=xfix ( {xlBoolean} ,  0)  ; 
sector3=xfix ( {xlBoolean} ,  0)  ; 
sector4=xfix ( {xlBoolean} , 0) ; 
sector5=xfix ( {xlBoolean} , 0) ; 
sector6=xfix ( {xlBoolean} ,  0)  ; 
z=xf ix ( { xlUnsigned, 10,0}, tempv) ; 
elseif  tempv<=2'^171-l 

sectorl=xfix ( {xlBoolean} ,  0)  ; 
sector2=xfix ( {xlBoolean} ,  1)  ; 
sector3=xfix ( {xlBoolean} , 0) ; 
sector4=xfix ( {xlBoolean} ,  0) ; 
sector5=xfix ( {xlBoolean} ,  0) ; 
sector6=xfix ( {xlBoolean} ,  0)  ; 
z=xf ix ( {xlUnsigned, 10,0}, tempv-171) ; 
elseif  tempv<=3'^171-l 

sectorl=xfix ( {xlBoolean} ,  0)  ; 
sector2=xfix ( {xlBoolean} ,  0)  ; 
sector3=xfix ( {xlBoolean} , 1) ; 
sector4=xfix ( {xlBoolean} ,  0) ; 
sector5=xfix ( {xlBoolean} ,  0) ; 
sector6=xfix ( {xlBoolean} ,  0)  ; 
z=xf ix  ( {xlUnsigned,  10,0},  tempv-2'^171 )  ; 
elseif  tempv<=4'*'171-l 

sectorl=xfix ( {xlBoolean} , 0) ; 
sector2=xfix ( {xlBoolean} ,  0) ; 
sector3=xfix ( {xlBoolean} ,  0) ; 
sector4=xfix ( {xlBoolean} ,  1)  ; 
sector5=xfix ( {xlBoolean} , 0) ; 
sector6=xfix ( {xlBoolean} , 0) ; 
z=xf ix  ( {xlUnsigned,  10,0},  tempv-3'^171 )  ; 
elseif  tempv<=5'^171-l 

sectorl=xfix ( {xlBoolean} ,  0)  ; 
sector2=xfix ( {xlBoolean} ,  0)  ; 
sector3=xfix ( {xlBoolean} , 0) ; 
sector4=xfix ( {xlBoolean} ,  0) ; 
sector5=xfix ( {xlBoolean} ,  1) ; 
sector6=xfix ( {xlBoolean} ,  0)  ; 
z=xf ix  ( {xlUnsigned,  10,0},  tempv-4'^171 )  ; 

else 

sectorl=xfix ( {xlBoolean} ,  0)  ; 


sector6,  z] 
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end 


sector2=xfix ( {xlBoolean} ,  0)  ; 
sector3=xfix ( {xlBoolean} ,  0)  ; 
sector4=xfix ( {xlBoolean} ,  0) ; 
sector5=xfix ( {xlBoolean} ,  0) ; 
sector6=xfix ( {xlBoolean} ,  1) ; 
z=xf lx  ( {xlUn signed,  10,0},  tempv-S'^lVl )  ; 


ramp2mod.m 


function  z  =  ramp2 (x) 

gain=xfix ( {xl Signed, 20, 19}, 1/1800) 

z=xf  lx  ( { xl  Signed,  14 , 13  } ,  x'^gain)  ; 


thetaconv2.m 

function  [y]  =  thetaconv(x) 
gainl  =  xf  lx  ( {xlSigned,  14 , 10  } ,  2'^3 . 14 )  ; 
gain2  =  xf lx ( { xlSigned, 14 , 10 } , 1/gainl ) 
if  x<0 

y=xfix  ( {xlUn signed,  10,0},  (x+gainl)  '^gain2'^1024) 
else 

y=xfix  ( {xlUn signed,  10,0},x'^gain2'^1024)  ; 
end 
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c. 


MATLAB  M-FILES  FOR  CHIPSCOPE  INTERFACE 


Black_box_dc_machine2_config.m 

function  code_conf ig ( this_block) 

%  Revision  History: 

o, 

o 

%  18-Dec-2008  (15:15  hours): 

%  Original  code  was  machine  generated  by  Xilinx's  System 

Generator 

%  after  parsing  H:\Docs\work_files\Docs\classes\EC3130_2009\DC 

%  machine  lab\black_box_dc_machine2 . vhd 

this_block . setTopLevelLanguage ( ' VHDL ' ) ; 

this_block . setEntityName ( ’ code ' ) ; 

%  System  Generator  has  to  assume  that  your  entity  has  a 
combinational 
%  feed  through; 

%  if  it  doesn't,  then  comment  out  the  following  line: 
this_block . tagAs Combinational ; 

this_block . addSimulinkInport ( ' ind ' ) ; 
this_block . addSimulinkInport ( ' ila_clock ' ) ; 
this_block . addSimulinkInport ( ' ind2 '  )  ; 

this_block . addSimulinkOutport ( ' outd ' ) ; 
this_block . addSimulinkOutport ( ' open_loop ' ) ; 
this_block . addSimulinkOutport ( ' speed ' ) ; 

outd_port  =  this_block . port (' outd ') ; 
outd_port . setType ( ' UFix_l_0 ' ) ; 

open_loop_port  =  this_block . port ( ' open_loop ' ) ; 
open_loop_port . setType ( ' UFix_l_0 ' ) ; 
speed_port  =  this_block . port (' speed ') ; 
speed_port . setType ( ' UFix_8_0 ' ) ; 

o_ 

o 

if  ( this_block . inputTypesKnown) 

%  do  input  type  checking,  dynamic  output  type  and  generic  setup  in 
%  this  code  block. 

if  ( this_block . port (' ind '). width  1); 

this_block .  setError  (' Input  data  type  for  port  "ind"  must  have 
width=l . ' ) ; 
end 

this_block . port ( ' ind' ) . useHDLVector (false) ; 
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if  (this_block . port (' ila_clock '). width  1); 

this_block . setError (' Input  data  type  for  port  "ila_clock"  must 
have  width=l . ' ) ; 
end 

this_block . port ( ' ila_clock' ) . useHDLVector (false) ; 

if  (this_block . port (' ind2 '). width  48); 

this_block .  setError  (' Input  data  type  for  port  "ind2"  must  have 
width=4  8  .  '  )  ; 
end 

end  %  if  ( inputTypesKnown) 

o_ 

o 


if  ( this_block . inputRatesKnown) 
setup_as_single_rate ( this_block, 
end  %  if (inputRatesKnown) 

o, 

o 


elk' 


ce'  ) 


%  Add  addtional  source  files  as  needed. 

%  I -  I  Add  files  in  the  order  in  which  they  should  be 

%  compiled.  |  If  two  files  "a.vhd"  and  "b.vhd"  contain  the  entities 
%  entity_a  and  entity_b,  and  entity_a  contains  a  |  component  of  type 
%  entity_b,  the  correct  sequence  of  |  addFileO  calls  would  be:  | 

%  this  block . addFile (' b . vhd ') ;  |  this  block . addFile ('a.vhd'); 


%  this_block . addFile ( ' ' ) ;  this_block . addFile ( ' ' ) ; 

this_block . addFile ( ' black_box_dc_machine2 . vhd ' ) ; 

return ; 


o. 

o 

function  setup_as_single_rate (block, clkname, cename) 
inputRates  =  block . inputRates ; 
uniqueInputRates  =  unique ( inputRates ) ; 

if  ( length (uniqueInputRates ) ==1  &  uniqueInputRates ( 1 ) ==Inf) 

block . setError (' The  inputs  to  this  block  cannot  all  be  constant.'); 
return; 
end 

if  (uniqueInputRates (end)  ==  Inf) 
hasConstantInput  =  true; 

uniqueInputRates  =  uniqueInputRates ( 1 : end-1 ) ; 

end 

if  ( length (uniqueInputRates )  1) 


59 


block . setError (' The  inputs  to  this  block  must  run  at 
rate .  '  )  ; 

return; 

end 

theinputRate  =  uniqueinputRates ( 1 ) ; 
for  i  =  1 : block . numSimulinkOutports 

block. outport (1) . setRate (theinputRate) ; 

end 

block . addClkCEPair (clkname, cename, theinputRate) ; 
return; 


single 
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APPENDIX  C :  SIMULINK/  XILINX  MODEL  OF  WIND  TURBINE 
GENERATOR  SYSTEM  (DFIG-SIDE) 
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A/D  CONVERTER  (VOLTAGE) 
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Pogtive  (Pk)  and  Negative  (Nk)  bits  are  created 
based  on  Aktrangtionsand  Bk 
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rotor 

V.aUloVjdl  contraller 


run 


ROTOR  CONTROLLER 
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PIDs 


labc  to  iqde 
SPEED  PI 
CURRENT  Pis 
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APPENDIX  D:  SIMULINK/  XILINX  MODEL  OF  WIND  TURBINE 
GENERATOR  SYSTEM  (SUPPLY-SIDE) 
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Vabc  s THETA  CALCULATION 
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APPENDIX  E:  TRANSFORMATIONS 


^ab  ^bc  -^b~^c  -^b~ {~^a  ~^b)’  whcrc  =  0  for  a  wye  connected  capacitor  bank  (with  the  neutral  floating) 


Note  also  that  i^  +  4  +4^0  ^  wye  connected  transformer  winding  set. 


73 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


74 


LIST  OF  REFERENCES 


[1]  R.  Pena,  J.  C.  Clare,  and  G.  M.  Asher,  “Doubly  fed  induction  generator  using 
back-to-back  PWM  converters  and  its  application  to  variable  speed  wind-energy 
generation,”  lEE  Proc.-Electr.  Power  Appl.,  Vol.143,  No  3,  May  1996. 

[2]  Shuhui  Li  and  Timothy  A.  Haskew,  “Analysis  of  Decoupled  d-q  Vector  Control 
in  DFIG  Back-to-Back  PWM  Converter,”  Power  Engineering  Society  General 
Meeting,  2007,  IEEE,  24-28  June  2007,  pp.1-7. 

[3]  G.A.  Smith  and  K.A.  Nigim,  “Wind-energy  recovery  by  a  static  Scherbins 
induction  generator,”  lEE  Proc.,  Vol.128,  Pt.  C,  No  6,  November  1981. 

[4]  Jason  G.  Massey,  “Douby  Fed  Induction  Machine  Control  For  Wind  Energy 
Conservation,”  Master’s  Thesis,  Naval  Postgraduate  School,  2009. 

[5]  H.N.  Robey,  H.O.  Stevens  and  K.T.  Page,  “Application  of  Variable  Speed 
Constant  Frequency  Generators  to  Propulsion-Derived  Ship  Service,”  Naval 
Engineers  Journal,  May  1985;  pp.296-305. 

[6]  Joseph  E.  O'Conner,  "Field  Programmable  Gate  Array  Control  of  Power  Systems 
in  Graduate  Student  Laboratories,"  Master’s  Thesis,  Naval  Postgraduate  School, 
2008. 

[7]  Paul  C.  Krause,  Oleg  Wasynczuk,  and  Scott  D.  Sudhoff.  “Analysis  of  Electric 
Machinery  and  Drive  Systems,”  John  Wiley  &  Sons,  2002. 

[8]  Alexander  Julian,  Giovanna  Oriti,  “FPGA  Based  Digital  Implementation  of 
Naturally  Sampled  Space  Vector  Modulation,”  ECCE  2009,  to  be  presented,  20- 
24  Sep  2009. 

[9]  Alexander.  L.  Julian,  Notes  for  EC4130  (Advanced  Electrical  machinery 
Systems),  Naval  Postgraduate  School,  2007  (unpublished). 

[10]  Andrew  M.  La  Valley,  “Design  and  Implementation  of  a  Motor  Incremental  Shaft 
Encoder,”  Master’s  Thesis,  Naval  Postgraduate  School,  2008. 


75 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


16 


INITIAL  DISTRIBUTION  LIST 


1 .  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Dr.  Jeffrey  Knorr 

Electrical  Engineering  and  Computer  Department 
Code  EC/Ko 

Naval  Postgraduate  School 
Monterey,  California 

4.  Dr.  Alexander  Julian 

Electrical  Engineering  and  Computer  Department 
Code  EC/Jl 

Naval  Postgraduate  School 
Monterey,  California 

5.  Dr.  Roberto  Cristi 

Electrical  Engineering  and  Computer  Department 
Code  EC/Cl 

Naval  Postgraduate  School 
Monterey,  California 

6.  Capt.  Lynn  J.  Petersen  USN 
Oxford,  Pennsylvania 

7.  Steven  B.  Swindler,  PE 

Naval  Surface  Warfare  Center 
Carderock  Division  Code  9820 
West  Bethesda,  Maryland 

8.  Joe  Borricinni 
NAVSEA  PHL  Code  9820 
Philadelphia,  Pennsylvania 


77 


